Low- Temperature Magnetic Penetration Depth in rf-Wave Superconductors: 
Zero-Energy Bound State and Impurity Effects 
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. We report a theoretical study on the deviations of the Meissner penetration depth A(T) from its 

£\j ■ London value in d-wave superconductors at low temperatures. The difference arises from low-energy 

surface Andreev bound states. The temperature dependent penetration depth is shown to go through 
O , a minimum at the temperature T m o ~ y^o/\oT c if the broadening of the bound states is small. 

The minimum will straighten out when the broadening reaches T m o- The impurity scattering sets 
up the low-temperature anomalies of the penetration depth and destroys them when the mean free 
path is not sufficiently large. A phase transition to a state with a spontaneous surface supercurrent 
is investigated and its critical temperature determined in the absence of a subdominant channel 
activated at low temperatures near the surface. Nonlinear corrections from Andreev low-energy 
bound states to the penetration length are obtained and shown, on account of their broadening, to 
O ' be small in the Meissner state of strong type II superconductors. 

o 

Oh 1 I. INTRODUCTION 

The low-temperature behavior of the magnetic penetration length in d-wave superconductors is in general a great 
3 ■ deal more complicated than that of their isotropic s-wave cousins. The changing sign of the order parameter, ac- 
cording to where one looks on the Fermi surface, entails coherent zero-energV|-jwJftw-Fenergy bound states in d-wave 
superconductors localized at smooth or almost smooth surfaces or interfacesElaLrta 9>H. These bound states feature 
pecuUajj- 1 low-temp| Crati J | rq iC pntri . bu . tkm s t p the magnetic penetration lengthl3(see alsoQ) and the zero-bias conductance 



C i _ 

' A minimum in the penetration depth of Y B0L2CU3O7-S filmsS and grain boundary junctions^ was thus interpreted 
as evidence for low-energy Andreev bound states. A conventional shielding-current contribution to the Meissner effect 
■ would obviously just monotonically reduce the penetration depth when the temperature goes down. On the other 
hand, a paramagnetic contribution from low-energy bound states increases the penetration depth. The interplay of 
these two effects amounts to a minimum in the penetration depth as a function of the temperature. The characteristic 
j— ( ' temperature T m Q of this anomaly is shown to be the order yj (£o/Ao)T c <C T c if the broadening 7 of the bound states 
is sufficiently small. At this temperature region the bound state contribution to the penetration depth competes with 
lO ' the low-temperature correction from shielding supercurrents to its zero-temperature value. 

An alternative explanation of an upturn in the penetration depth is possible in compounds whose bulk para- 
magnetic properties grow when the temperature goes down, like in the j-eJ^cttpn-doped cuprate superconductor 
Nd\.^Ceo.ibCuOi-y. There the paramagnetism arises from Nd 3+ ional3'c2lEilE3. We will not discuss these com- 
pounds below. 

There is yet another important temperature associated with the magnetic penetration depth, T s ~ (£,o/^o)T c . If 
a given crystal orientation does not carry quasiparticle Andreev bound, states, a nonlocal effect can take over as a 
correction to the zero temperature penetration depth in the clean limited. Then in other orientations which do admit 
Andreev states, the bound-state contribution and the spontaneous surface supercurrent in particular, can-in turn 
overwhelm the nonlocal effect. At T ~ T s the bound state paramagnetic contribution to A in the clean limited is the 
order of the total London penetration depth Ao from the screening currents. In the absence of sub-dominant pairing 
chpHiels, a spontaneous surface supercurrent brought about by the bound states may arise below the temperature 
r s Eal23 (see alsoEZlon a similar effect of spontaneous magnetization brought about by low energy interface bound states). 
Having in mind high-temperature superconducting compounds, we will discuss strong type II superconductors. Then 
(£o/-V)) is easily the order 0.01, and the low temperature range splits up into at least three areas staked out by 0, T s 
and T m o (T s ~ (£o/Ao)T c <C T m o ~ \J (£o/Ao)T c <C T c ). Quasiparticle scattering off impurities or surface roughness 
and inelastic processes may also play an important role if they bring about a broadening 7 of the bound states the 
order or greater than the characteristic temperature T m Q (T s ). 

We assume below that nonmagnetic impurities dominate the scattering and the broadening. Nonmagnetic impurities 
in superconductors with an anisotropic order parameter are known to be pair breaking. They suppress T c analogously 
to what happens to isotropic superconductors with magnetic impurities. Assuming superconductors always clean 
within the conventional definition £ <C I, we disregard this kind of effects throughout the article. Even then impurity 
broadening of Andreev bound states in anisotropically paired superconductors can be significant. Since the broadening 
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removes singularities in the density of states (for instance, 5-peaks from quasiparticle bound states) as well asin other 
related physical quantities, superconductors can be sensitive to extremely small concentrations of impuritiesEj. This 
is analogous to the role of pair breaking and small anisotropy of the gap in the Riedel anomaly in isotropic s-wave 
superconductors. The Riedel anomaly is associated with the BCS singularity in the density of states. Pair breaking 
and small anisotropy of the gap are known to wipe out the BCS singularity in the density of states averaged over the 
Fermi surface, and control the height of the Riedel peakc3. 

The emphasis of the present work is on the various effects of broadening on the low temperature anomalies of the 
Meissner effect. The zero-energy pole-like term of what is known as the quasiclassical Green's function was exploited 
in the investigation. Broadening is introduced into the pole-like term simply sliding the pole along the imaginary 
energy axis. With small broadening, relatively simple expressions are found for the penetration length in the two 
lowest-temperature regions defined above. If T s / m Q\ ^5 7 *C T c , the growing 7 can wipe out the low-temperature 
anomalies. Beginning with the critical broadening y s ( m ), anomalies at T g ( m ) are fully destroyed. It turns out that 
unitary scatterers need to come with significantly larger scattering rates r s ( m ) than Born impurities in order to achieve 
the critical broadening 7 s ( m ). This effect is peculiar of the impact of impurities on the Andreev bound states as seen 
in the local density of states and Josephson critical currentsc3. For this reason, the requirements the mean free path 
must meet for the low temperature anomalies to show up are sensitive to the strength of the impurity potential and 
very different in the unitary and the Born limits. 

For Born scatterers, the shortest normal-state impurity mean free path I which preserves the low temperature 
upturn at T ~ T m is shown to be Ao ^$ I. This looks quite restrictive althoug h conc e ivabl y compatible with the 
strikingly large low temperature mean free paths in some high-T c compoundsl^r 1 lrnrTP 4 rP 5 l. For the spontaneous 
surface supercurrent in the absence of a subdominant component at the surface, we find the threshold Ag/£o < I. This 
demands extraordinary clean samples not available for the time being. On the other hand, the requirements set by 
unitary scatterers are much weaker and probably can be met. In this case surface roughness is likely to control the 
broadening and the experimental observability of the effects. 

We also examine what the Andreev bound states do to the nonlinear Meissner effect. AtJow temperatures T <C T c , 
the field Hq at which the nonlinear response of the bound states saturates in the clean limitcil is much weaker than the 
one from the screening current. Ignoring the broadening, Hq is a linear function of the temperature. With T < T s , 
nonlinear effects from the bound states become important already in the Meissner state. Close to the transition to 
the state with a spontaneous surface supercurrent, a nonlinear term entering into the Landau mean-field free energy 
is important also in a weak external field. The broadening 7 introduces another field, H 7 characterizing the nonlinear 
consequences of the bound states at irT < 7. For sufficient broadening irT s <C 7, we get H c \ <C H 1 and the nonlinear 
terms are shown always to be small in the Meissner state. 



II. THE UPTURN IN THE LOW-TEMPERATURE DEPENDENCE OF THE PENETRATION DEPTH 

Our considerations are based on the quasiclassical matrix Green function which describes quasiparticle excitations 
in thermal equilibrium. The quasiclassical propagator g(pj, x, e n ) satisfies Eilenberger's equations, which have a 2 x 2 
particle-hole matrix form 

f ie n + -v/-A(R)j f 3 - A(p/,R) - a(pf,R;e n ),g(p f ,TL;e n ) 

-Mv/-Vr <?(p/,R;e n ) = , (1) 
g 2 { Pf ,R;e n ) = -ir 2 i , (2) 

where e n = (2n + \)ttT are the Matsubara energies, p/ the momentum on the Fermi surface, v/ the Fermi velocity, 
A the vector potential, A the order parameter matrix, and a the impurity self-energy. A symbol with a 'hat' denotes 
a matrix in the Nambu space. 

The propagator g and the order parameter matrix A parameterize as 

The gap function A(p/,R) is related to the anomalous Green function / and must be determined self-consistently. 
The diagonal part g(p/,R, e n ) of the full matrix propagator g carries information on the electrical current density 

j(R) = 2eTN f J2(v f g(Pf,Re n )) Sf . (4) 
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Here Nf is the normal state density of states per spin direction and < ... >g, means averaging over quasiparticle 
states at the Fermi surface. 

Let an anisotropic singlet strong type II superconductor occupy the right half-space x > 0. A magnetic field 
is applied along the z-axis. The induced supercurrent and the vector potential (in the gauge divA(R) = and 
vanishing in the bulk) have only y-components. The linear current-field relation in general has a nonlocal form, i.e. 
j(x) = - f+°°Q(x,x',T)A(x')dx'. 

For strongly type II superconductors with nodes.i±i the order parameter, a nonlocal current-field relation can be of 
importance only at very low temperatures T < Tj?3. Hence, a study of the penetration depth at low temperatures 
T s <C T ~ T m o <C T c may be carried out disregarding nonlocal effects. Then a magnetic field enters into Eq.(^|) only 
together with the Matsubara frequencies e n — i^Vf jV A(x) in the argument of the Green's function. The kernel Q(x, T) 
can then be written 

«,,n-2^2£ £ ( v y P;) ^£te)) s . (5) 

n— — oo 

In the presence of zero-energy surface bound states, the pole-like term in the propagator becomes dominating at 
temperatures r « T c . Surface bound states as well as their paramagnetic response are localized on the scale of the 
coherence length at the surface, however, while the conventional screening current has an avenue of the huge thickness 
of the penetration depth. That is why the zero-energy bound-state contribution to the penetration depth remains a 
small low-temperature correction to Ao = Xq(T = 0) at all temperatures T>r s (in particular, at T ~ T m o). The 
contribution from surface bound states must be viewed together with a low temperature correction from the screening 
current as small low-temperature imports to the zero temperature London penetration depth Ao- Then the total 
kernel of the form Q(x,T) = c . 2 + SQ(x,T) includes only the lowest order corrections in 5Q(x,T). 

. o 

Solving the Maxwell equation 

1 s 4tt 

perturbatively with respect to the last term delivers a first order approximation to the vector potential 



A"{x) - -^A(x) SQ(x, T)A(x) = (6) 

An c 



A(x)=A^(0) 



-\-oo 

'<>'i-£)-^/^(-^H'M-£ 



(7) 



The kernel SQ(x 1 T) incorporates only a contribution from the bound states and a low-temperature correction from 
the screening current. 

The penetration depth is defined as A = H{x)dx/ H{Q) = — A(0)/A'(0). Expanding this to first order in 5Q 
and extracting the low temperature correction from the screening current for the case of a superconductor with a line 
of nodes 

A(T) = A + aA ^ - ^ / Q bound (x, T)dx . (8) 
T c c J 
o 

Here a is a coefficient of the order of unity which depends on the shape of the Fermi surface and on an angular slope of 
the order parameter near the nodes. For instance, for a quasi-two-dimensional d x 2_ y 2 tetragonal superconductor with 
a cylindrical Fermi surface (with a principal axis z) and order parameter A(0) = Ao cos(2</> — 2a), one gets a w 0.32. 

Kernel Q bound (x,T) takes negative values. It is a paramagnetic contribution from zero-energy bound states to 
Eq.(|). One obtains Q bound (x,T) from Eq.(||) substituting instead of the full expression for g(pf,x,s n ) only its 
singular part (pole-like term) g s (pf,x,e n ). Associated with zero energy surface bound states, this term vanishes in 
the bulk on the scale of the coherence length £o- It has longer tails only towards the nodes. Node contributions do 
not dominate, however, in the following expressions. The presence of zero-energy surface bound states is crucial in 
the reasoning. All sectors of the Fermi surface associated with a sign change of the order parameter in a quasiparticle 
reflection from the surface, contribute significantly to the results. This allows us to neglect, to a good accuracy, the 
factor exp(— 2x/Ao) under the integral sign in Eq.(|§|). 

The analytic expression for the pole-like term has been found in the clean limit and for a smooth surface in Rcf. 
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g.<p f ,* t e n ) = — ^ q) + ^ Q)| e( P/ ) exp [-^^ J q \*<Pf,*)\&) • (9) 
The effective surface order parameter |A(p/,0)| introduced in Eq.(^J), is denned 

cxp I - | — rr / |A(p/,a;')|da;' I da; . (10) 



|A( P/ ,0)| \vf, x (Pf)\J \ WfAPf)\ 
o \ o 



Here we distinguish between incoming pj and outgoing quasiparticle momenta in a reflection event. For specular 

reflection, the momentum parallel to the interface is conserved. Function 0(f>/) is equal to unity where zero energy 
bound states occur on the Fermi surface (i.e. where the order parameter in the bulk taken for incoming pf and 
outgoing p momentum directions have opposite signs), and vanishes elsewhere. 

Substituting Eq.fl) in Eq.(|^), one can easily sum over the Matsubara frequencies. Integration over the space 
coordinate x in Eq.flq) then yields the penetration depth: 

KT) = Ao + a^-X Q + ^f^ 20 (vl y (p f )\v Lx (Pf)\Q(Pf)) Sf , T c ^ « T « T c . (11) 

For a three dimensional superconductor with a spherical Fermi surface one has the relation A 2 , = 3c 2 /(87re 2 v 2 iVf ). 

Then the coefficient in front of the third term in Eq.dnl) is ^ 7T \> . Analogously, for a simple model of a quasi-two- 

8T V/ 

dimensional superconductor with a cylindrical Fermi surface, A 2 , = c 2 /(47re 2 v 2 iV/) and the coefficient ^ $ . 
In particular, for a d x 2_ y 2-wave superconductor with a cylindrical Fermi surface, we get from Eq.(|i"T|) 

A(T) = A + a— A + ^||sin 3 /?|-| COS 3 /3|| , T c -^- «T«T C , (12) 
T c 6T 1 1 A 

where f3 = a + (7r/4) is the angle between the surface normal and the direction to a node of the order parameter, 
while a is the angle between the surface normal and the crystalline a-axis along its positive lobe. 



We note that the correction from zero energy bound states to the penetration depth (the third term in Eq.(ll)) 
has a quite universal form. It is independent both of the spatial profile of the order parameter near a surface and 
its particular anisotropic structure (basis functions). Therefore, this correction depends only on the type of pairing, 
which determines regions on the Fermi surface with opposite signs of the order parameter. For example, expression 
( |l2| ) is valid irrespective of a particular form of a momentum direction dependence of the basis function for a d-wave 
order parameter of given symmetry. 

The ratio of a supercurrent density at the surface j h ° und (x — 0, T) to the one j s (x^ r , T) at a characteristic distance 
x scr (£o ^ x scr Ao) from the surface can be estimated for a clean superconductorcil at T <C T c and a smooth surface 
as \j b ° und {x = 0,T)/j s (x scr ,T)\ - 4n\l\Q bound (x = 0,T)|/c~ T c /T. This verifies that at low temperatures T«T C 
the paramagnetic current j h ° und {x,T) dominates over the shielding current near the surface within a relatively small 
characteristic scale £o- 

The temperature dependent terms in Eq.(|l"l"|) behave in very different fashions from each other. They come from the 
conventional shielding currents and from the zero-energy bound states. Growing with decreasing temperature, the dia- 
magnetic screening currents monotonically reduce the penetration depth. On the other hand, Andreev surface-bound 
states respond paramagnetically and increase the penetration depth when the temperature goes down. Disregarding 
the broadening effects, Eq.(jll|) delivers the following estimate for the field of the low-temperature minimum of the 
penetration depth: 

(13) 

where £ is of the order of unity for crystalline orientations with sufficient amount of momentum directions admitting 

1 /2 

zero-energy bound states. Otherwise £ is a small quantity. For a d-wave superconductor Q oc ||sin 3 /3| — |cos 3 /3|| 
and vanishes for (3 = 45° (i.e. for a = 0), when there are no zero energy bound states. 

Broadening of the bound states can substantially modify the conditions for the presence of a minimum in the low 
temperature dependence of the penetration depth. For a small broadening j(pf) T c we simply replace the factor 
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■jJ- in the expression Eq.(||) for the pole- like term with t — — — j^-t -, — r 

Eq. ( |ll| ) is generalized to the following form: 



Taking into account the broadening 



A(T) = A + 



a-A 



2e 2 N f X 2 



c 2 T 



v f,y 



( P f)WfAPf)MPfW 



i(pf) \ 

2ttT ) 



(14) 



Sf 



Here and below ip( x ) is the digamma function and ip'i x ) ~ its derivative. 

Eq. (jl4|) is a reasonable representation of the role of a broadening in the low temperature anomaly of the penetration 
depth. The minimum lies at T m o « 1.8y£o/AoT c for momentum independent broadening in a ^2 „ J( 2 -superconductor 
in the clean limit 7 <C 7rT with the orientation a — 45°. With increasing broadening it drifts to lower temperatures 
(becoming less pronounced at the same time) till T mi w 0.4^/£ /AoT c at 7 w 0.96T m o, where it evaporates. As an 
example, the low-temperature correction to the penetration depth is shown in Fig. 1 in the vicinity of T mj for various 
values of the momentum independent broadening. 
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FIG. 1. Low-temperature correction to the penetration depth (in units of Ao) in a d, x ?_ y i superconductor with a cylindrical 
Fermi surface and the orientation a — 45°. The temperature is measured in units of T c . The parameter fo/Ao is chosen to be 
w 0.01, where £0 = ( v //27rTc). The curves are given for three values of the broadening: 7 = 0.10T C (dashed line), 7 = 0.15T C 
(solid line) and 7 = 0.19T C (dashed-dotted line). 



There are various contributions to the broadening of the bound states associated, in particular, with surface rough- 
ness, nonmagnetic and magnetic impurities and inelastic scattering. We now pin-point the origin of the broadening, 
assuming that nonmagnetic impurities dominate the scene. 

IT 

With Born scatterers ji, y -=£ (see Ref. |2^) and the coefficient of the order of unity can be estimated within the 
simple model of spatially constant order parameter. Then we easily get the shortest normal-state impurity mean free 
path I which admits a low-temperature upturn: Ao ^$ I- In high-temperature superconductors one should distinguish 
between I and the actual mean free path in the normal state at T c incorporating significant contributions from inelastic 
processes. Impurity scattering dominates there at low temperatures already in the superconducting state where the 
collapse of inelastic scattering takes place. For instance, below 20 K in \lS6to£ij* 3 .pj'_£ there is a regime of extremely 
long and weakly temperature dependent quasiparticle scattering timed 3 "Ir tPtP^ 3 tE3 usually interpreted as due to 
feeble impurity scattering in high- purity samples. 

For scatterers with sufficient strength of impurity potential there are practically no restrictions on the impurity 
scattering rate in contrast to what was found above for Born impurities. For unitary scatterers with scattering rates 
r„ <C T c the broadening of the zero-energy bound states is exponentially smaller: j u = By/A T U exp (— 6A /r„). 
A scattering rate r u which leads to a given broadening 7„ is almost independent of a constant coefficient B in the 
pre-exponential factor, while it is sensitive to the model dependent parameter b in the argument of the exponential 
function. Within the simple model considered in Ref. |2^, one gets b ~ 1. 

For temperatures T J$ y/T u Ao, the share of the penetration depth from the shielding currents must be modified for 
unitary scatterers. This leads instead of the linear term in Eqs.(||), (|ll|), (jlj), ( flij ) to a quadratic low temperature 
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correction of the form XqT 2 / (r 1 / 2 A^ 2 ^J to within a factor the order of unityliEE. Correspondingly, T m Q given 
in Eq.(|l3|) is valid for unitary scatterers only if T m Q > ^/T u Aq, which sets an upper limit on the scattering rate: 

The T 2 -term instead of the linear one in Eq.fO) delivers an estimate for the location T m d of the low temperature 
minimum of the penetration depth modified by unitary scatterers: 

T md » f|^V /3 (T u Al) 1/e , T u > f^A . (15) 
\Ao/ A 

This expression replaces Eq.([l3|) in the case of unitary scatterers with the scattering rate T u > Aq. The minimum 

slowly drifts to higher temperatures with increasing T u . It does not melt away at any T u <C T c . The normal-state 
impurity mean free path must just be large on the scale of the coherence length. 

We conclude that observation of the low temperature upturn of the penetration depth in samples with I < Ao 
is evidence for both Andreev bound states and a sufficiently large strength of the bulk impurity potential in the 
superconducting compounds. For unitary impurities one needs to take into account the broadening that arises from 
surface roughness which then very probably controls the total broadening. The same effect with Born scatterers 
demands the normal-state impurity mean free path larger than the London penetration depth. 



III. ZERO ENERGY BOUND STATES AND SPONTANEOUS SURFACE CURRENT 



Throughout this section the broadening of the zero-energy bound states is assumed small. We look at a clean d-wave 
superconductor with a smooth surface. Its crystal-to-surface orientation shall admit zero-energy surface bound states 
and feature an upturn in the penetration depth. Below the upturn temperature T m o, imagine a great deal of space 
for A to grow, firstly as described by the perturbative result Eq. (j]~T|) . Then a secontLpOrder phase transition occurs 
at T ~ T s -C T m Q into a state which carries a spontaneous surface supercurrentE§n3E3. We shall find an analytic 
expression for the transition temperature and discuss the impact of impurities on the effect. The transition implies 
the absence of subdominant channels activated at low temperatures close to the surface on account of the presumably 
large surface pair breaking in, tlie||fiLpminant component of the order parameter. Otherwise a spontaneous current can 
arise at higher temperaturesuSclnj. 

There is experimental evidenced] for a phase transition on the (110) surface in Y Ba2Cu^OT-s at T = IK ^> 
(£o/Ao)Tc-n It was interpreted as associated with an activated near surface subdominant channel of the order 
parameter!! Eot-some other crystal-to-surface orientations, however, a subdominant component can be not present 
near a surfaceQ'ta. Zero energy bound states can still arise for a noticeable part of quasiparticle trajectories. Our 
theoretical study is relevant to these cases. 

In order to find an equation for the transition temperature, one has to admit a paramagnetic contribution to the 
penetration depth at least as large as the diamagnetic one. Then a perturbation treatment of the preceding section 
is not adequate. In this context we develop an approach based on the integral form of Eq.([|) and take into account 
only the terms in SQlx, T) brought about by the bound states. In other words, a contribution only from the pole-like 
term Eq.@ needs to be taken into account in Eqs.(||) for the kernel which enters into Eq.(||). The kernel SQ(x,T) 
varies on the characteristic scale £o and is associated in the clean limited with large contributions to the magnetic 
field at the surface at temperatures T < (£o/Ao)T c . We therefore disregard the nonlocal temperature correction from 
the Meissner current to SQ(x,T). 

We transform Eq.(^|) into the integral form 



A{x) = 



A(0) 



2vrAn 



+ 00 



dx'Q bound (x',T)A(x') e A o 



x 

Ao. 



-Hoc- 



2 "' V ' dx'Q bound (x',T)A(x') 



x — X 



(16) 



The two terms on the right hand side of this equation obey very different scales. The first decays exponentially 
in the depth on the scale Ao while the last term vanishes for x 3> £0 along with the kernel Q bound (x,T). The 
kernel Q bound (x — 0,T) can be estimated (see preceding section) for a clean superconductor and a smooth surface 
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as 2n\lQ bound (x = 0,T)/c - T c /T. Then, in accordance with Eq.(|l6|), the approximate formula ( 1 - 



A{0) * 

^qT c /AqT is established for a relative deviation of the vector potential A(x scr ) taken at the distance x scr (£ <C x scr <C 
Ao) from its value A(0) at the surface. The deviation reflecting the bound state contribution to the vector potential 
turns out to be small at all temperatures T ^S> (£oMo) ^c> m particular, for T ~ T s ~ (£o/Ao)T c . Varying on the 
scale small terms in the expression for the vector potential at temperatures T ~ T s are of importance only when 
differentiating A(x). After that they can already noticeably contribute to the expression for the magnetic field. 
Indeed, a spatial differentiation of Eq.(|l6|) leads to 



H(x) = 



A(0)-^^ / dx'Q fcotind (x',T)A(a;') ( e^o 



+00 , x , 



x 



c 



dx'Q 6oun,, (« , ,r)^(!C / ) U A ° +e A ° • (17) 



The second term in the square brackets remains neglig ibly small - (^T c /AgT) A(0) < A(0) as compared with A(0) 
for T > (Co/Ao) T c . The last term of Eq.@ is the order of (£ T C /A T) (A(0)/A ). For a deviation of the magnetic 

f H(x ) \ 

field at x = x scr from its value at x — 0: — 07 nT — 1 I ~ (£oTc/AqT). Hence, the bound state contribution to 



, #(0) ^ 

the magnetic field can be viewed as a small perturbation as compared with the shielding contribution unless T < T s . 
Considering (£oMo) TV^- ~ we can discard the second term in the square brackets but have to keep track of 
the last term in Eq . (|l7|) . Choosing x = in Eq.([T^), the small terms x'/\q in the exponential functions under the 
integral sign can be taken to vanish. For the same reason and within the same accuracy one can treat the vector 
potential under the integral sign in Eq.(|l7|) as constant in space A(0) discarding small terms in the vector potential 
which vary on the scale £o- All this results in an explicit relation between A(0) and H(0) and therefore 

A= . (18) 



l + lZ^£y Q b ° md (x,T)dx 


Proceeding like in the derivation of Eq.([ll]) above, we find that the paramagnetic (negative) sign of Q bound leads 
to a divergence of A at the temperature 



Ts = ^ (vL(P/)|v/,x(p/)|e( P/ )) S/ . (19) 



c 



For the model c?-wave superconductor with a cylindrical Fermi surface one gets from Eq.(p^) 

3A 



r s -^||sin 3 /3|-|c OS 3 /3||T c , (20) 



where = v//27rT c . 

The divergence of A implies the existence of a nontrivial solution to Eq.([l6|) in a vanishing external magnetic 
field. Indeed, if we let H(0) = 0, A(0) 7^ 0, then Eq.(p7|) transforms, with the same approximation as above, into 
the relation: 1 = — (AtiXq/c) / + °° Q bound (x', T)dx' , which defines the transition temperature T s into a state with a 
spontaneous surface supercurrent. 

The nontrivial solution at T s is a result of interplay between the paramagnetic supercurrent which originates in the 
zero energy bound states localized within £0 on the one hand and the diamagnetic supercurrent distributed over the 

region x ~ Ao on the other. The latter compensates for the magnetic field from the bound states at the surface in order 

+00 

to satisfy the boundary conditions in the absence of an external magnetic field. Then J j(x)dx = always applies 



being a consequence of the full screening of the spontaneous surface magnetic field in the bulk of a superconductor. 

Under these conditions the Bloch theorem, in general, admits spontaneous surface currents!^. The magnetic part of 

+00 



a superconducting free energy J 



A' 2 {x) + ^Q{x,T)A 2 (x) 



dx vanishes at T — T s and becomes negative below 
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T s on account of negative sign of the paramagnetic kernel Q bound (Gibbs and Helmholtz free energies coincide in zero 
external magnetic field). The result is an energetically favorable state with a spontaneous surface supercurrent below 
Ts. 

The broadening 7 of the bound states modifies Eq.(|19|): 



2e 2 N f X 



VfJPf)\vfAPf)\®(PfW 



1 , jiPfY ] 

2 2vrT s J 



(21) 



Sf 



The broadening prevents the appearance of a spontaneous surface current unless 7 < J^~T C - This is a very strong 
restriction. If Born impurities control the broadening, they admit spontaneous surface supercurrent only with ex- 
tremely large values of the mean free path Aq/£o ~ lOOAo ^ I, unrealistic for high temperature superconductors. 



Unitary scatterers impose a much weaker restriction T u < 26Aq / In 



A 



£oMV6>) 



0.1T r . Then, however, surface 



roughness probably dominates the broadening and can destroy the state with a spontaneous surface supercurrent. 



IV. NONLINEAR MEISSNER EFFECT FROM LOW ENERGY BOUND STATES 

It is important in the derivation of Eq. ([l8|) that the kernel Q bound varies much faster in space than the screening 
currents. Then contributions of the paramagnetic current carried by surface Andreev-bound states at temperatures 
(£o/Ao) 2 T" c <C T, can result in significant spatial variations of the magnetic field near the surface while in the weakly 
spatially dependent vector potential. This leads to Eq.([l8j) on the basis of the local current-field relation. 

It is straightforward to show within the same framework that a nonlinear penetration depth X n i (T, H) incorporating 
contributions both from screening currents and from zero-energy bound states is described as 

X nl (T,H) = Ka(T,Hscr) ^ (22) 



4ir\ s n T(T,H scr ) 



J Q b n °r d (^T)dx 



where X^f (T, H scr ) is a contribution from screening supercurrents to X n i(T,H), taken at an effective value of the 

. +00 +00 / * 

field H scr = H{0) - J H{x)dx J Q b ° und (x' ,T)dx' . Here H(0) is the external magnetic field. The second term 
^ ' 

describes the field of the zero-energy bound states inside the superconductor at distances x ~ x scr (£0 x scr <C Ao), 
as can be seen in Eq. (Jl7|) . A paramagnetic response of zero-energy bound states {Q h °i und < 0) increases the field 
to be screened by diamagnetic supercurrents (H scr = H(x scr ) > H(fi)). This leads, in general, to more pronounced 
nonlinear terms in A^f* (T, H scr ) as compared to disregarding the contribution from zero-energy bound states. In the 
case of spontaneous surface supercurrent H scr differs from zero even in the absence of an external field. We assume 
the condition H scr < H c i for the Meissner state to be stable in the magnetic field on account of a paramagnetic 
influence of the bound states. 

Small nonlocal low temperature corrections to the penetration depth from screening currents can be taken into 
account in Eq.([2^) as perturbations to A^ r (T, H). For a nonlocal current-field relation a penetration depth X n i(T, H) 
is actually a functional of the spatial profile of the magnetic field. 

Nonlinear corrections from the shielding supercurrent to the Meissner effect can be given in terms of the dimen- 
sionless ratio p — (H/Ho), where Ho is usually the order of the thermodynamic critical field ~ &o/(Xo£o). Hence, 
they are always small in strong type II superconductors in the Meissner state. In isotropic s-wave superconductors, 
the first nonlinear correction to the penetration depth cx p 2 . In superconductors with nodes in the order parameter 
(for instance, d-wave ) a term linear in p can arise for particular crystal orientations at low temperaturesCJO. The 



linear term, however, is quite sensitive, tp nonlocal effectsEa and the impurity influence, in particluar, at sufficiently 
large strength of impurity potentialsOCJ. 

A nonlinearity in the magnetic response of low energy Andreev surface bound states has, in general, a very different 
field scale Hq. In a clean limit Hq(T) = (<&oT/Xvf), where A is determined by Eq. (p2[) . At low temperatures T«T C 



one always gets H„(T) < H . For instance, H (T m0 ) - ^ &H - 0.1iJ , H (T S ) - ^fH cl < H cl - O.Olffo- 

Moreover, at sufficiently low temperatures T <^ T s , i.e. well below the transition to the state with a spontaneous 
surface supercurrent, a paramagnetic response from the bound states may become seriously nonlinear already in the 
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Meissner state Hq(T) ~ H < H c ic3/Bl. We will show, however, that the broadening of the bound states introduces a 
new held scale H-y — y^-^o coming into play at 7rT < 7. For 7 ^> T c - 0.01T C nonlinear corrections from Andrcev 

low-energy bound states to the penetration length turn out always to be small in the Meissner state, even at T = 0. 

As a pole-like term Eq.(^) decays exponentially on the scale ~ £0 for almost all momentum directions admitting 
bound states, we consider a local nonlinear current-held relation 

j{x) = 2eTN f J2(vf9 s (Pf,x,£n-i^yf,yA(x))) (23) 

S„ f 

for the current via the bound states. One can also set the vector potential A(x) in the kernel equal to A(x = 0). 
Then we easily generalize the reasoning in the derivation of the third term in Eq. (O) . Substituting into Eq.(p3|) the 
expression Eq.(H) for the pole- like term with the pole shifted in accordance with the broadening, we find: 

h , ieNf I (l 7(P/)+*;v / ,,(p / )A(0)\ 
Q% und (x,T)dx = -rd(v f , y (p f )\v f , x (Pf)MPfW o + ^ ■ ( 24 ) 



In Eq.(|2^) the Fermi surface is assumed symmetric in reflections across the xz-plane. Then averages over the Fermi 
surface of odd powers of Vf iV vanish, no matter whether they are multiplied by \vf,x(Pf)\®(Pf) or not. This applies, 
in particular, to a d-wave superconductor with a cylindrical Fermi surface whose principal axis z is parallel to the 
boundary for arbitrary orientations of the two other crystal axes xo, yo- 

If ^VfA(0) -C max(27rT, 7), one can expand the -0-function in Eq.(|24|) in powers of the small parameter 
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AA^f and bound states AA* ; to be small, one can represent them in the first approximation as additive contributions 
to the total penetration depth \ n i(T, H) w A(T) + AA^ r + AA^ ; . The nonlinear correction from screening currents 



min -ff(O) / Hq(T), H(0) j H 1 I . Considering nonlinear corrections to the penetration depth from screening currents 



takes the form AA*<f = 



I ' ^^nl 



Quantities A(T) and Aq(T) being the zero-held val- 



ues of X n i(T,H), A^ r (T, H) respectively, satisfy Eq.(nq). Bound states renormalize nonlinear response from screening 
currents already in this approximation. Thus, the explicit analysis of AA^ r can be done combining the results of the 
preceding section and Refs. [y], [l2], Apart from too close to the transition temperature T s , the nonlinear correction 
to the penetration depth from the bound states is: 



AA* ; 



In the limit ^VfA(0) ^> 2-kT when the argument of the ^-function in Eq. (p4|) is large, we obtain 

00 

\^T)dx = -^L (v f , y (p f )W f , x (Pf)\e(Pf)™ctan i ^^^ )) ■ ( 26 ) 



-\bound f 
Inl 



Then the broadening rather than the temperature fixes the bound state contribution to the penetration depth. 
As shown above, there is no state with a spontaneous surface current with 7 ^> T s . Then A ~ Ao and evfA(0) j ay <C 

H(p)/H c ±. Since H(0) < H c i in the Meissner state we estimate evfA(0) j cy <C 1 and obtain in this limit from Eq.(|26|) 

AA , = 4^a| ( v ^ (p/) | v ^ (p/) |e( P/ )) S/ - 

3cV ( y fJPf)\ v fAPf)\®(Pf)) Sf ■ (27) 

For a given A^ r (T, H ) , Eq. (]2^) in general should be solved with respect to A„; (T, H ) in accordance with Eq. (|4|) , 
since A(0) = — A n ;(T, H)H(0). This is particularly important close to the transition temperature T g , where the Landau 
theory of second order phase transitions is applicable. Then hrst nonlinear term turns out to be the order of the 
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zero-field paramagnetic contribution in the denominator in Eq.(|22j). Ignoring a weak field dependence of X^f(T,H) 
stipulated by screening currents, we obtain from Eqs.(^2|), (|24|) the following equation for X n i(T,H): 

(Jr - lj «A„;(T, H) + r/iJ 2 A^ ; (T, H) = X s n f(T s ) , (28) 

where H = H(Q), 



T=T 3 



T s is described by Eq. (||l]) and T s (T) is the result of the substitution T s — > T on the right hand side of Eq. ( |2~l| ) . The 
broadening is assumed to be sufficiently small for admitting the phase transition. 

The role of the order parameter in the phase transition can be played by a surface magnetization ms = 

+00 +00 +00 

/ (M(x) — Moo) dx = j; J dxxj s (x) = -J- J H(x)dx — -J-\ nl (T, H)H (0), which, for simplicity, we choose con- 
77 7r 

stant in space along a smooth surface. The magnetization M enters by the conventional definition j = c • curlM, and 
Afoo = —H(0)/4tt. Then the Landau free energy per unit surface Ts which leads to the same equation for ms as 
implied in Eq.(^8|) has the form 

Fs = & ■ Qr - l) rn% + pm% - m s H , (29) 

where a — 2nK,/\ s ^ l r (T s ), (3 — 16Tr 3 rj/X^ r (T s ), H is the external field. As for a conventional order parameter in a 
strong field near T s one gets ms(T s , H) oc iJ 1 / 3 , which entails X n i{T s , H) oc H~ 2 / 3 . 

Finally, in the limit of very small broadening 7 -C ^^-T c , Eqs.(|22]) and ( p6| ) give 



1^(0)1 



AA^ - X nl (T,H) - KT(T,H) = ' u (\vf,y(pf)vfAPf)\®(Pf))s f ( 30 ) 



at temperatures ■^■T' c «T« —h-^-T c . Since AA^ is at least the order of Ao, we put here A*f (T, iJ) w Ao disregarding 

small nonlinear corrections from the screening currents in the Meissner state. The approximate inverse proportionality 
of the penetration length the magnetic field implies the presence of a spontaneous surface magnetization weakly 
dependent on H . 



V. CONCLUSION 



We have examined the paramagnetic contribution from surface zero-energy Andreev bound states to the low- 
temperature penetration length of d- wave superconductors in the Meissner state. The paramagnetic current is localized 
within several coherence lengths near the surface and grows larger in the clean limit when the temperature goes down. 
A broadening of the bound states chokes their contribution and determines their actual role in shaping the penetration 
length. We found that the upturn in the low temperature penetration depth lies at T m Q ~ y/£o/ XqT c in the clean 
limit where the paramagnetic contribution from the bound states can be handled with perturbation theory same as 
small low-temperature corrections to the penetration depth from the screening current. The minimum broadening 
capable of straightening out the upturn is 7 w T mo . 

Furthermore, we examined the penetration depth when the bound states must be kept track of beyond perturbation 
theory. A divergence of A(T) was found at the phase transition to a state with spontaneous surface supercurrent. This 
transition occurs only with smallish broadening, 7 < (£o/Ao)T c . In the clean limitEa and at low temperatures, there 
is a nonlinear regime of the paramagnetic current already in magnetic fields substantially weaker than the fields for 
the nonlinear effects to show up in response of shielding supercurrents. The broadening of the bound states modifies 
and weakens the nonlinear response. 

Specifying an origin of the broadening as associated with nonmagnetic impurity scattering, we obtained restrictions 
on the mean free path admitting the low temperature anomalies. The conditions turn out to be sensitive to the 
strength of the impurity potential and very different in the unitary and in the Born limits. The Born impurities are 
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shown to easily prevent the anomalies of the penetration depth taking place at least well below T m o. By contrast, 
unitary scatterers with sufficiently small normal-state scattering rate T u <C T c admit the transition to a state with 
spontaneous surface supercurrent at T s ~ (£q/\q)T c . In the latter case, however, surface roughness very probably 
dominates the broadening and controls the bound state contribution to the low-temperature penetration length. 
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